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Abstract 

The forward-backward multiplicity correlation strength is calculated for 
arbitrary nucleus-nucleus collision in the framework of the wounded nucleon 
model. Discussion of our results in the context of the recent STAR data in 
AuAu collisions at ^/s = 200 GeV is presented. It is suggested that the ob- 
served (i) growth of the correlation coefficient with centrality and (ii) approx- 
imately flat pseudorapidity dependence of the correlation strength for central 
collisions are due to the fluctuations of the number of wounded nucleons at a 
given centrality bin. 

PACS: 25.75.-q, 25.75.Gz 
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1 Introduction 



Recently the STAR collaboration announced the results [1] on the forward-backward 
multiplicity correlations in nucleus-nucleus collisions. It was found that the corre- 
lation strength (defined below) was larger than in an elementary proton-proton 
collisions and it remains constant (at least for the most central collisions) across 
the measured midrapidity region. This result was interpreted in the framework of 
the color glass condensate [2] or dual parton [3] models, which suggests the possible 
formation of high density partonic matter in central AuAu collisions at ^/s = 200 
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GeV. Other theoretical investigations concerning the problem of forward-backward 
multiplicity correlations in hadronic collisions can be found in Refs. [4-12]. 

The main difficulty, however, is to distinguish between correlations arising from 
the presence of the quark-gluon plasma and correlations that do not depend on 
this new phenomenon. These need to be understood, controlled and subtracted 
in order to access the true signal of the quark-gluon plasma. The natural ground 
to study this problem is the wounded nucleon model [13]. Indeed, it is the sim- 
plest superposition model in which a nucleus-nucleus collision is constructed from 
an elementary nucleon-nucleon collisions. More precisely, the number of produced 
particles in nucleus-nucleus collision is proportional to the number of wounded nu- 
cleons, i.e., nucleons that underwent at least one inelastic collision. A Monte Carlo 
analysis of this problem in the very simplified wounded nucleon model was already 
presented in Ref. [4]. 

Our main conclusion is that the STAR data [1] can be naturally understood in the 
wounded nucleon model (at least for the most central collisions). We conclude that 
the observed growth of the correlation coeflicient with centrality and approximately 
fiat pseudorapidity dependence of the correlation strength are due to the fiuctuations 
of the number of wounded nucleons at a given centrality bin. 

The correlation coefficient (or correlation strength) b is defined as 



where ub and np are event by event particle multiplicities in backward B and for- 
ward F pseudorapidity^ intervals, respectively. The main ingredients which allow to 
evaluate b in the wounded nucleon model are (i) recently obtained pseudorapidity 
particle density from a wounded nucleon p{r)) and (ii) particle multiplicity distri- 
butions measured in pp collisions in different forward and backward intervals. The 
fragmentation function p{rj), shown in Fig. 1, was obtained by analysing the PHO- 
BOS data on dAu collisions [14] at ^/s = 200 GeV in the wounded nucleon model 
[15] and the wounded quark-diquark modeP [16]. In a completely independent way, 
mainly based on the recent NA49 collaboration data [17], analogous wounded nu- 
cleon fragmentation function was constructed in Ref. [18]. For the multipficity 
distributions measured in pp collisions we take the negative binomial (NB) fits [19] 



where n is the average multiplicity and 1/k measures deviation from Poisson distri- 
bution. 

^Our discussion is valid for any longitudinal variable, not necessarily pseudorapidity. 

^In this case p (77) = 1.2^(77) + 0.8U{i]), where F{r]) and U{r]) are the particle densities from 
wounded and unwounded constituents, respectively. 
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Figure 1: A wounded nucleon fragmentation function at = 200 GeV cm. energy. 

In the next section the correlation coefficient b for the collision of two arbitrary 
nuclei is derived. In section 3 we focus on the collision of two symmetric nuclei and 
look closer at the midrapidity and fragmentation regions, where b can be written in 
a particularly simple form. Our results are discussed in the context of the recent 
STAR data in section 4 and section 5 where also some comments are included. In 
the last section our conclusions are listed. 



2 Model 

The problem is to calculate the correlation coefficient (1) in two given pseudorapidity 
regions B and F under assumption that the contribution to the multiplicity in these 
two intervals is provided by independent contributions from left- and right-moving 
wounded nucleons. It is similar to the assumption of independent hadronization of 
strings in the dual parton model [3]. The picture of independent left- and right- 
moving sources of particles is the main assumption of the wounded nucleon model. 
There are many phenomenological and experimental evidences supporting this idea 
[15,18,20]. 

It is convenient to construct the generating function 

H {ZB, zf)=Y,P K, np) zl^z''/, (3) 

where P {nBinp) is the probability to find ub particles in B and in F. In general 
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we may have many sources of particles (wounded nucleons), thus P {ub, rip) can be 
expressed as 



P{nB,nF)^ ^ W {wL,WR)P{nB,nF;wL,WR) , 



(4) 



WL,WR 



where W {wl, wr) is the probabihty distribution of the numbers of wounded nucleons 
moving left wl and right wr, respectively. P {nB.np\WL^WR) is the probability to 
find hb particles in B and n^? in F under condition of and wr wounded nucleons 
in left- and right-moving nucleus, respectively. 

As derived in the Appendix, the generating function (3) reads 

-kwL/2 



H{zb,zf)= ^ W {wl,wr) ^1 + ^\pLBil - zb) + Plf {1 - Zf)]^ 



X 



X 



|l + ^ \Prb (1 - zb) +Prf (1 - zf)]^ 



-kwji/2 



(5) 



where prf denotes the probability that a particle originating from the right-moving 
wounded nucleon goes to F interval, under the condition that this particle was found 
either in _B or F (and analogous for prb,Plb and plf)- These probabilities satisfy 
natural conditions 

Plb + Plf = 1, Pre + Prf = 1- (6) 

These numbers can be easily calculated. Indeed, they depend only on positions and 
sizes of B and F as well as the shape of the wounded nucleon fragmentation function 
p (r/) . For instance, prf has the form 

JpPiv) drj 



Prf = 



(7) 



L+Fpi^) dV 

The parameters n and k come from the NB distribution fit (2) to the pp multiplicity 
distribution data in the combined interval B + F. These parameters are well known 
for various energies and different pseudorapidity intervals [19]. Moreover 



n = 2 



P iv) dr}. 



(8) 



B+F 



It is worth to notice that formula (5) contains all information about the multiplicities 
in B and F, as well as their dependence on the number of wounded nucleons. 

Using definitions (1) and (3) we obtain [6 = U/D\: 

d^H dH dH' 



U 



D 



OzbOzf 



d^H dH 
+ 



dzB Ozf 
dH 



dzl 



dzF \dzF 



2" 



(9) 
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Performing appropriate differentiations we obtain 



4C/ 

^ = PlbPlf 



+ PrbPrf 



+ {plbPrf + PlfPrb) [{wlWr) - (wl) (wr)] , 



+ 

(10) 



and 



4D 



Plf 



/ 2\ / \2 , 2(wl) 



/ 2 \ / \2 . 2 {wr) 
{Wr)-{wr) +— ^ 



+ 



+ '^PlfPrf [{wlWr) - (wl) (wr)] + 2 



Plf {wl) + Prf (wr) 



n 



(11) 



In the above expressions (...) represents the average with respect to W {wl,wr). 
For instance (wl) is the average number of wounded nucleons in the left- moving 
nucleus. 



3 Fully symmetric case 

The result presented in the previous section is valid for any colliding nuclei. In case 
of symmetric collisions we of course have (wl) = (wr) and (wl) = (w'j^). Moreover, 
studying correlations in symmetric (around 7] = 0) intervals, i.e., piB = Prf = P 
and Plf = Pre = 1 — p, we obtain 

^ _ [K) - {wRf + 2 {wr) /k] + C2 [{wlWr) - {wRf] 

C2 [{wl) - {wr)' + 2 {wr) /k] + Ci [{wlWr) - {wr)'] + 2 {wr) In' 

where 

Ci = 2p(l-p), C2-I-C1. (13) 

It is worth to notice that formula (12) simplifies for two cases. 

(i) Midrapidity. Considering two narrow rapidity intervals B and F around 77 = 
we have p ~ 0.5. It leads to a particularly simple expression 

-1 



6 = 1 



n/2 (■»'-)-(»')' 



1 + v v + 



4 I A; {w) 



(14) 



where w = ^ wr is the number of wounded nucleons in both colliding nuclei. 
This formula allows to notice the growth of b with increasing scaled variance of the 
number of wounded nucleons (Jw — {w)]'") / {w) . 

(ii) Fragmentation region. Assuming that intervals B and F are separated 
enough so that F can be populated only by right-moving wounded nucleons and 
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B only by the left- moving ones, that is p = 1, we obtain 

^ ^ {wlWr) - {wr) 

In this case b > only due to the fluctuations of the number of wounded nucleons, 
i.e., 6 = if (wlWr) = {wr}^. 

This closes the theoretical discussion of the problem. 



2 2 
n k 



— _L 



(15) 



4 Results 

Recently the STAR collaboration presented results [1] on correlation coefficient b 
for AuAu collisions at = 200 GeV. The backward B = (-^ - 0.1, + 0.1) 
and forward F = — 0.1,-^ + 0.1) intervals of width 0.2 each were located 
symmetrically around r] = with the distance At] between bin centres ranging from 
0.2 to 1.8 with an interval of 0.2. The measurement was performed for different 
centrality classes defined via the number of produced particles in the central rapidity 
region.'^ 

As argued in Ref. [22] different centrality selections (e.g., via impact parameter, 
number of wounded nucleons, number of produced particles) give the same average 
number of wounded nucleons (w). However, as was shown in Ref. [4], they lead to 
rather different Q, = [{w'^) — {w)'^]/ (w), except the most central collisions, where Q 
weakly depends on the centrality class definition. In consequence, direct comparison 
of our result (12) with the STAR data can be performed only for the most central 
collisions. In case of non-central collisions the comparison is not straightforward. 
Indeed, wounded nucleon model does not describe correctly the multiplicities in 
AuAu collisions, thus not allowing to impose experimental centrality class cuts on 
the number of produced particles. 

We performed our calculations with the centrality class definition via the number 
of wounded nucleons w — wl + wr in both colliding nuclei (obviously the impact 
parameter fluctuations are also included). We performed Monte-Carlo calculations 
[23] for five centrality class selections: - 10%, 10 - 20%, 20 - 30%, 30 - 40% 
and 40 - 50% what correspond to w > 275, 275 > w > 197, 197 > w > 139, 
139 > w > 94 and 94 > w > 60, respectively. The corresponding results for (wr), 
(wj^), (wlWr) and Q are presented in Table 1. In our MC calculations for the nuclear 
density profile we took the standard Woods-Saxon approximation with the nuclear 
radius R = 6.38 fm and the skin depth d = 0.535 fm [24]. For the nucleon-nucleon 

^For instance — 10% centrality class corresponds to events with the number of produced 
particles (in the central region) larger then 430 [21]. 
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% 




{<) 


(wlWr) 


n 


0-10 


163 


26841 


26793 


3.04 


10-20 


116.8 


13803 


13723 


2.07 


20-30 


83.1 


7015 


6938 


1.71 


30-40 


57.5 


3378 


3315 


1.40 


40-50 


37.83 


1478 


1432 


1.26 



Table 1: The MC results for (wr), (wj^), {tVLWRj and Q = [{w^) — (w)'^]/ (w) 
corresponding to five 10% centrality classes defined via the number of wounded 
nucleons w = Wi, + wr in both colliding nuclei. 

interaction profile we used the black disk approximation"^, i.e., the interaction takes 
place only if the transverse distance between two coUiding nucleons is smaller than 
y^cr/TT, with the total inelastic pp cross section cr = 42 mb. 

In Fig. 2 the calculated correlation coefficient h (12) for — 10% and 10 — 20% vs. 
the distance A77 between bin centres is compared with the STAR data [21]. Taking 
Eq. (7) into account we obtain p = 0.51, 0.52, 0.55, 0.56, 0.58 for Ar^ = 0.2, 0.6, 1.0, 
1.4, 1.8, respectively. NB distribution fits to pp multiplicity data in the midrapidity 
region give approximately constant n = 0.96 (central plateau) and k = 1.8 [19].^ It is 
interesting to note that for the — 10% most central events, where direct comparison 
with the data is possible, the wounded nucleon model can explain more than 85% 
of the effect. 

In Fig. 3 the correlation coefficient b (12) for 20 - 30%, 30 - 40% and 40 - 50% 
centrality events vs. the distance Arj between bin centres is shown. The wounded 
nucleon model predicts larger values of b than observed, however, as explained at 
the beginning of this section in case of non-central collisions the direct comparison 
with the data cannot be performed [unknown precise value of [(tf^) — (w)"^]/ (w), see 
Eq. (14)]. The main experimental finding, however, that the correlation coefficient 
is approximately flat^ in the midrapidity region as a function of At] is very well 
reproduced in the wounded nucleon model. This feature of the model can be easily 
understood from Eq. (14). Indeed, n and k are approximately constant (central 
plateau) in the midrapidity region and the value of p is close to 0.5, which is a 
consequence of the longitudinal structure of the wounded nucleon fragmentation 

^We also performed calculations for the Gaussian approximation. We observe the week depen- 
dence of our results on the pp interaction profile. 

''From Ref. [19] it may be concluded that k is slightly increasing to fc « 2 for Arj = 1.8. This 
effect, however, practically does not influence numerical values of b. 

^This is not true for the 40 — 50% centrality events. We will come back to this point in the next 
section. 
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Figure 2: The STAR data points compared with the results of the wounded nucleon 
model for the correlation coefficient b for two most central events vs. the distance 
Arj between bin centres. The width of each bin equals 0.2. 

function, shown in Fig. 1. 

In Fig. 4 the correlation coefficient b (12) for - 10%, 10 - 20% and 30 - 40% 
centrality events vs. the distance A77 in the broader range is shown. Taking Eqs. 
(7) and (8) we obtain {p, n) = (0.7, 0.97), (0.85, 0.79), (1, 0.53) for Ar] = 4, 6, 
8, respectively. The values of parameter k are not known precisely within these 
intervals, however, as can be concluded from [19] they should not be larger then 
k — A. As shown in Fig. 4, where the results are presented for k — 1.8 and 
A; = 4 for Arj > 4, this uncertainty practically does not influence our flnal results. 
The reduction of the correlation coefficient b at A77 = 8 (at this point 6 for 
peripheral collisions) is fully determined by the suppression of particle production 
from a wounded nucleon to the backward hemisphere^, see Fig. 1. 



5 Comments 

Following comments are in order. 

(a) It is well-know that the wounded nucleon model signiflcantly underestimates 
the multiplicities in AuAu collisions [25] . Contrary to the model assumption multi- 

''For instance, assuming that the contribution from a wounded nucleon is symmetric around 
T] = (i.e. p = 0.5 at any Ar]) we would obtain b « 0.2 at Arj = 8 for 30 — 40% centrality events. 
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Figure 3: The STAR data and the results of the wounded nucleon model for the 
correlation coefficient b for non-central events vs. the distance A77 between bin 
centres. As explained in the text in this case direct comparison between data and 
model cannot be performed. 



plicity from a wounded nucleon depends on the number of collisions it underwent. 
In order to take this effect into account wc multiplied n and k by the ratio 7 of 
the measured multiplicity in AuAu collisions [25] to the prediction of the wounded 
nucleon model nw/2 [13]. For the — 10% most central collisions we approximately 
obtain 7 ^ 1.6 and consequently b ~ 0.59, which is in very good agreement with the 
measured value, see Fig. 2. 

(b) As seen in Figs. 2, 3 and from Eq. (14) the correlation coefficient calculated 
in the wounded nucleon model is always flat in the midrapidity region in contrast to 
the 30 — 40%, 40 ^ 50% peripheral AuAu or pp collisions. In the present paper we 
suggest that the fluctuation of the number of wounded nucleons may be responsible 
for the large value of the forward-backward correlation coefficient in central AuAu 
collisions. For peripheral collisions, however, this source of correlations is becoming 
less important (the value of [(w^) — (w)^]/ (w) decreases) and obviously the mech- 
anism responsible for correlations in elementary pp collisions play a major role. In 
our approach we cannot describe the precise shape of the correlation coefficient in pp 
collisions^ as a function of A77, thus at the same time our approach is not applicable 
to the peripheral AuAu collisions. 

^Neglecting fluctuations in Eq. (14) we obtain in the midrapidity region the constant value of 
0.2. 
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Figure 4: Wounded nucleon model prediction for the correlation coefficient b in the 
broad range of the distance A77 between bin centres. The width of each bin equals 

0. 2. 

(c) Encouraged by the success of our approach we also provide the prediction 
for the correlation coefficient b for the most central PbPb collisions at the LHC 
energy ^/s = 5500 GeV. Performing appropriate MC calculations^ (described in the 
previous section) we obtained for the — 10% centrahty events the following value 
of [{w'^) — (w)"^]/ (w) — 2.55. Once the parameters n and k are measured at a given 
B + F interval in pp collisions, the predictions for the correlation coefficient b in 
the midrapidity region can be easily obtained from Eq. (14). As an example we 
present the result for b in the midrapidity region with the forward and backward 
intervals identical to those at the STAR measurement. The needed parameters 
n ~ 1.7 [26] and k ~ 1.5 [19,27] arc taken from extrapolations to the LHC energy. 
In consequence we obtain b ~ 0.6 or 6 pa 0.7 if the correction to the wounded nucleon 
model, discussed at the beginning of this section, is taken into account. Our value 
is close to the prediction reported in Ref. [6]. However, in our approach the large 
value of b is only due to the fluctuation in the number of wounded nucleons, which 
if neglected (it corresponds to pp coUisions) we obtain b ^ 0.35 in contrast to the 
value reported in Ref. [6]. 

(d) The wounded nucleon pseudorapidity fragmentation function, shown in Fig. 

1, extends far beyond its own hemisphere. As discussed in the previous section this 
feature is partially responsible for the approximately constant value of b in the midra- 

^Here R = 6.62 fm, d = 0.546 fm and = 67 mb [24]. 
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pidity region. It is interesting to note that similar longitudinal structure is present 
in the dual parton model (DPM) [3] , where the long longitudinally extended strings 
are stretched between quarks and diquarks of the projectile and target, respectively. 
In general, models that can explain the long-range forward-backward correlations 
are models that introduce long extended objects in rapidity [2, 3, 5-7, 10]. Moreover, 
in DPM the growth of the correlation coefficient is due to the fluctuations in the 
number of elementary inelastic collisions, which is similar to the fluctuations of the 
number of wounded nucleons present in our approach. Therefore it is not surprising 
that the two models lead to similar qualitative results [1]. However, the wounded 
nucleon model is in better agreement with data. 

(e) Similar longitudinal structure provides the QCD inspired color glass conden- 
sate model (CGC) [28], which includes many features of DPM. In this approach [2] 
the long extended color fiux tubes and the fluctuations of the number of gluons allow 
to understand the main features of the STAR data. Moreover, it was shown recently 
[29] that the soft ridge structure observed at RHIC [30] can be naturally understood 
in the CGC/glasma motivated phenomenology, which is rather difficult to obtain in 
the framework of the wounded nucleon model. This problem is currently under our 
investigation. 

(f ) It would be interesting to perform similar calculation of the correlation coeffi- 
cient b in the framework of the wounded quark-diquark model [16, 31], which proved 
to be quite successful in description of particle production in pp, dAu, CuCu and 
AuAu collisions. In this model the number of produced particles is proportional 
to the number of wounded quarks and diquarks, which are assumed to be the con- 
stituents of each nucleon. Here the growth of the correlation coefficient is due to 
the fluctuations of the number of wounded quarks and diquarks at a given centrality 
bin. 

(g) In the present approach we implicitly assume that particles are produced 
directly from wounded nucleons. It would be interesting to check an effect of in- 
termediate resonances (clusters) production. We expect this effect to influence the 
forward-backward multiplicity correlations in the midrapidity region for peripheral 
AuAu and pp collisions. 

6 Conclusions 

Our conclusions can be formulated as follows. 

(i) We have studied the forward-backward multiplicity correlations in the frame- 
work of the wounded nucleon model [13]. In this model particles are produced 
independently from the left- and right-moving nucleons that interacted in inelastic 
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way at least once. An analytical expression for the correlation coefficient (strength) 
for the collision of two arbitrary nuclei and at any forward F and backward B 
intervals was derived. 

(ii) The main ingredients of our approach are: recently obtained long extended in 
pseudorapidity wounded nucleon fragmentation function [15, 16] and the multiplicity 
distributions measured in proton-proton collisions described by a negative binomial 
distribution. 

(iii) In the midrapidity region correlation coefficient can be written in a particu- 
larly simple form (14). This expression allows to explain the growth of the correla- 
tion coefficient with increasing scaled variance ([w — (tf)]^) / (w) of the number of 
wounded nucleons w in both colliding nuclei. 

(iv) We have performed explicit calculations for AuAu coUisions at ^/s — 200 
GeV. The backward/forward intervals were chosen according to the recent STAR 
measurement. Growth of the correlation coefficient with centrality as well as almost 
no pseudorapidity dependence in the midrapidity region was observed. Our results 
are in good qualitative agreement with the STAR data, although exact comparison 
can be performed only for the most central collisions. 

(v) Finally, predictions for the values of the correlation coefficient in the broad 
range of pseudorapidity were presented. 
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A Appendix: Generating function 

Let Phi'TiLB-i nip) be the probability that a left-moving wounded nucleon contributes 
ulb particles into B and nip particles into F interval [and analogous distribution 
PR{nRB,TiRp) for a right-moving source]. The probability to find Ub = n^B + i^rb 
particles in B and np = n^p + tibf particles in F in case of one left- and one 
right-moving wounded nucleon is given by 

PiuB^np)^ ^ PL{nLB,nLF)PR{nRB,nRp)5ll^^^^Jll^^^^^, (16) 

nLB,'nLF 
nRB,nRF 
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and the generating function (1,1 means one left- and one right-moving wounded 
nucleon) 

H {ZB, ZF] 1,1)= J] P (riB, np) zl^ z""/ = Hl {zb, zf) Hr {zb, zf) , (17) 

with 



Hl (zb, Zf) = ^ PiiriLB, tilf)^ 

nLB,nLF 

Hr{zb,zf)= J2 PR{^RB.nRF)z7''z'i,^^. (18) 



It is easy to check that in case oi wl left- moving and wr right- moving wounded 
nucleons we obtain 

H {zB, ZF\ WL, Wr) = [Hl {zr, Zf)]""' [Hr {zr, Zf)]'''' . (19) 

Suppose that Pi (n) is the multiplicity distribution from a single wounded nucleon 
in the combined interval B + F. Then 

PLiuLB, ulf) = Pi{n = ULB + ^^^) ^^^^ + ^^^) ' (p^^)"-- (plf)""" , (20) 

ULB^-riLF^- 

where plb and plf are defined in section 2. In consequence 

Hl {zb, Zf) = V Pi(n) [plbZb + PlfZfT ■ (21) 
Performing analogous calculations for the right moving part 

Hr (zb, Zf) = y] Pi{n) [prrZb + PrfZfT ■ (22) 

' 'n 

Assuming that the multiplicity distribution measured in proton-proton collision 
is described by a NB distribution with n and k (in the combined interval B + F), it 
is easy to show that Pi{n) (from a wounded nucleon) is given by a NB distribution 
with n/2 and k/2 [32]. Then using 

E„^'.("r=(i+^)"'". (23) 

we obtain 

H {zb, Zf; wl, wr) = |i ^ [1 - plbZb - PlfZf] | 

X 1 1 + - [1 - prbZb - PrfZf] I ■ (24) 



-kwL/2 

X 

-kwR/2 



Summing over W{wl^ Wr), i.e., the probability distribution of the number of wounded 
nucleons, and taking (6) into account we finally obtain (5). 
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